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ABSTRACT 

Many models for coronal loops have difficulty explaining the observed EUV 
brightness of the transition region, which is often significantly less than theoretical 
models predict. This discrepancy has been addressed by a variety of approaches in- 
cluding filling factors and time-dependent heating, with varying degrees of success. 
Here we focus on an effect that has been ignored so far: the absorption of EUV light 
with wavelengths below 912 A by the resonance continua of neutral hydrogen and he- 
lium. Such absorption is expected to occur in the low-lying transition region of hot, 
active region loops, that is co-located with cool chromospheric features and called 
"moss" as a result of the reticulated appearance resulting from the absorption. We use 
co-temporal and co-spatial spectroheliograms obtained with SOHO/SUMER and Hin- 
ode/EIS of Fe XII 1242A, 195 A and 186.88 A, and compare the density determination 
from the 186/195 A line ratio to that resulting from the 195/1242 A line ratio. We find 
that while coronal loops have compatible density values from these two line pairs, 
upper transition region moss has conflicting density determinations. This discrepancy 
can be resolved by taking into account significant absorption of 195 A emission caused 
by the chromospheric inclusions in the moss. We find that the amount of absorption 
is generally of order a factor of 2. We compare to numerical models and show that 
the observed effect is well reproduced by 3D radiative MHD models of the transition 
region and corona. We use STEREO A/B data of the same active region and find that 
increased angles between line-of-sight and local vertical cause additional absorption. 
Our determination of the amount of chromospheric absorption of TR emission can be 
used to better constrain coronal heating models. 
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Introduction 



The nature of the transition region continues to be the subject of significant debate. As part 
of the interface between the photosphere and corona, it plays an important role in the transport 
and deposition of the non-thermal energy believed to power the heating of both chromosphere and 
corona. It is now well established that the interface between the photosphere and corona does not 
resemble the view suggested by hydrostatic, spherically sy mmetrical an d time-averaged models. 
Instead, the chromosp here is a highly dynamic environment (Ruttenll2006) with jets, such as fibrils. 



ties dj 
et all 



spicu l es and mottles (Hanst een et al J 120061 : iDe Pontieu et al. 
2007l : bePontieu 



2007 



a; 



Rouppe van der Voort et al. 



2007bl) continuously perturbing the thermodynamic quantities and density 



stratification from their equilibrium values. As a result of these significant deviations from hy- 
drostatic equilibrium, the transition region is molded by the same dynamic spicular features as 
the chromosphere. There is ample evidence that between a height range of roughly 2,000 and 
10,000 km above the photosphere, plasma with temperatures ranging from 10,000 K to 1 MK 
co-exist and "mingle" in a highly dynamic fashion. However, the availability of time-averaged, 
one-dimensional models of the chromosphere and transition region have allowed coronal loop 
mod els to ignore th e complexities and inherently dynamic nature of this interface (for a review, 
see lKlimchukll2006r) . Neglecting the intricacies of the chromosphere can be problematic since the 
chromospheric heating mechanism requires 50 times more energy than the coronal heating mech- 
anism. Consequently, it is perhaps not surprising that a full explana tion of coronal heating re mains 
outstanding, with many unresolved issues and problems (but see lAschwanden et al J 120071 for a 
different approach). 

One of these issues that has been problematic for theoretical models of coronal loops is the 
large discrepancy between the observed and predicted emission from the transition region in coro- 
nal loops. A variety of theoretical models predicts emission in transition region lines much in 
excess of what is observed with EUV imagers such as SOHO/EIT and TRACE. This manifests 
itself in large differences in predicted brightness between active region loops and transition re- 
gion moss. Moss is low-lying, strong EUV emission that was discovered i n TRACE EUV images 



in the Fe ix/x 171 A and Fe xn 195 A passba nds (IBerger et al.lll999bl : be Pontieu et al.lll999l : 
Fletcher & de Pontieu 19991: Berger et . al.l 1 1999ab. It is the upper t ransition region at the base of 



hot, high pressure loops feerger et al.l 1 1999bh iMartens et al.l 120001) . Observations show that in a 
typical Fe IX/X 171 A EUV image of an active region, coronal loops and moss regions are roughly 
of the same brightness. Despite a variety of approaches, all theoretical models continue to pre- 
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dict much higher EUV emission from moss regions than from coronal loops. For example, most 
coronal loop models based on a steady heating mecha nism predict moss emission that is one to sev- 
eral orders of magnitude larger than what is observed (|Schrijver et al.ll2004c IWarren & Winebarger 
20061 : Iwarren et alJl2008h. _ While the introd uction of models based on time-dependent nanoflare 
heating (e.g. IPatsourakos & Klimchukl 12008) on average shows a smaller discrepancy, they still 
predict moss emission that is 10 times brighter than what is sometimes observed. Other models 
remove the discrepancy by introducing a filling factor in the transition region as a free parameter 
Jwinebarger et alJl2008h . 



However, none of these models take into account the fact that the transition region emis- 
sion in moss regions as observed by EUV imagers such as TRACE and EIT undergoes significant 
absorption by neutral gas through Lyman continuum absorption. Radiation of EUV emission be- 
low the head of the hyd rogen Lyman continuum at 912 A can be absorbed by neutral hydrogen 
(lAnzer & Heinzell 120051) . In addition, neutral helium can absorb EUV radiation with wavelengths 
shorter than 504 A and singly ionized helium can contribute to absorption for wavelengths shorter 
than 228 A. The effects of neutral hydrogen and helium and singly i onized helium on EUV ra- 
diation were s t udied extensively in the past dSchmahl & Orrall Il979l: boschek & Feldmanll 19821 
Mariskal[l992l baw et all Il995h I Judge et all 1 19951 : knzer & Heinzell I2OO5I) for a variety of solar 
regions including quiet Sun, coronal holes, prominences, the limb, and active regions, often with 
confusing results and interpretations. Part of this confusion likely stems from the fact that the 
amount of neutral plasma in observations of such a wide variety of features is highly variable and 
most often not directly observed, but simply inferred from the EUV emission. This is not the 
case for moss: high resolution observations made in the EUV (with TRACE) and chromospheric 
lines such as Ha (with ground-base d telescopes) clearly show the presence of neutral plasma at 
heights where EUV emission occurs (|Berger et al.lll999b|) . In fact, moss shows up as a reticulated 
pattern above plage regions precisely because the EUV emission from the 1-1.5 MK upper tran- 
sition region plasma is absorbed by mostly neutral gas in dense, chromospheric jets or fibrils that 
occur at similar heights and in the immediate surroundings of the upper transition region plasma 
(be Pontieu et alJll999h . 



The mingling of chromospheric and TR plasma at low heights has been confirm ed by sev- 
eral observational studies in the past few years. For example. be Pontieu et al. ( 2003b|) found that 
in moss regions the linecenter intensity of one of the dominant chromospheric lines (Ha) is very 
well correlated, both spatially and temporally, with Fe ix/x EUV emission observed with TRACE. 
Such a correlation is difficult t o understand unless chrom ospheric and TR plasma occur at similar 
heights. This is confirmed by iDe Pontieu et al.l (|2003al) who found that oscillations in moss, ob- 
served in EUV with TRACE, are associated with oscillations in the wings of Ha. These correlated 
oscillations can be explained by absorption of EUV emission by chromospheric fibrils that quasi- 
periodically rise high enough so that they can obscure the lower-lying upper TR emission from 
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neighboring field lines (IDe Pontieu et al.ll2004|) . Such a scenario is supported by extensive obser- 



that drive chromospheric plasma upward to reach heights of several 1,000 km ( 


3e Pontieu et al. 


2004 


De Pontieu & Erdelvi 


2006; 


Hansteen et al. 


2006; 


De Pontieu et al. 


2007a 


Heaaland et al. 


2007 


). 



An accurate determination of this absorption in moss regions is thus clearly an important in- 
gredient in determining the accuracy (or constraining parameters such as filling factors) of coronal 
loop models. In this paper we focus on determining what fraction of the EUV emission in moss 
regions is absorbed, so that future work can directly compare predicted intensities with observed 
intensities that have been corrected for this absorption. 

We use the density sensitivity of emission lines (of the same ion) formed above and below 
the Lyman continuum and com pare the obtained densities fro m several line pairs. We describe 
the observations with Hinode (IKosugi et al.1 120071) and S OHO jP oming o et al.N 1995D and the co- 
alignment of the EIS dCulhane et alJl2007h and SUMER dwilhelm et all 1995b spectra in §2. In §3 
we describe the discrepancies in density measurements for coronal loops and moss regions, and 
estimate t he amount of absorp tion of TR emission. We also observe the same moss region with 
SECCHI dHoward et alJbood) onboard the STEREO A and B spacecraft (kaiser et alJl2008h and 
use the dependence on viewing angle of the TR emission to estimate the center-to-limb variation 
of the absorption (§4). We demonstrate that such absorption compares well with what we detect 
in numerical simulations (§5). We finish by summarizing our results and briefly discussing the 
potential impact of our results on coronal loop modelling in §6. 



2. Details of EIS/SUMER Observations 

To study the differential absorption of EUV emission with wavelengths below and above 912 
A, we focus on emission lines from the same ion that is formed, under ionization equilibrium 
conditions, at temperatures that are high enough to clearly show moss (logT = 5.7 — 6.3). These 
conditions are quite stringent since most high-temperature lines above 912 A are forbidden lines 
and quite weak. In addition, the small-scale nature of moss structuring requires 1-2" class spectra 
in both lines. The launch of Hinode and its EIS instrument has opened up access to a range of 
Fe XII lines (186.88 and 195.1 A in particular) that together with SUMER's Fe XII 1242 A line 
provide the tools we require to address the absorption of EUV moss emission by neutral hydrogen 
and helium (and singly ionized helium). We have chosen these lines because the 186.88 A and 
195. 1 A line pair is sensitive to densities in the range of log n e = 8 — 10.5, an ideal range for active 
regions. In addition, the ratio of 195 to 1242 A is also sensitive to density, so that we can directly 
compare the densities derived from both line pairs. 
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We use an EIS raster of a weak active region (which lacks a NOAA number) taken from 
16:44-17:50 UTC on 14-Nov-2007 (see panels A, B of Fig. [[])• The active region was observed 
at a heliocentric position of (+150", -100"). The raster consists of spectra of a range of EUV 
emission lines (including Fe XII 186.88 and 195.1 A) at 256 positions (spaced 1" apart in solar x) 
observed from west to east at a cadence of 15 s. The data is despiked (i.e., cosmic ray removal) and 
calibrated using eis_prep . pro in the EIS solarsoft software tree. For each location in the raster, 
we calculate the total intensity in the Fe XII lines by summing over a range of spectral pixels that 
includes the full line, and subtracting an estimate of the continuum emission in a spectral region 
in the immediat e spectral vicinity th at is free of emission lines. The absolute calibration of EIS 



is described in iCulhane et al.1 (|2007r) and is estimated to be accurate to about 20%. We remove 
the pointing offset between the EIS rasters by using a cross-correlation technique between the two 
Fe XII rasters. 

During the same time frame SUMER observed this active region in a range of UV lines 
that include the Fe XII 1242 A emission line. Similarly to both EIS lines, it has a peak in the 
contribution function at log T = 6.15. The raster was taken from east to west and took much longer 
to complete (from 17: 12 UTC through 20:35 UTC) because of the different SUMER sensitivity and 
the weakness of the Fe XII line. The raster steps are 1.125" with three exposures of 30 s at each 
raster position to ensure an adequate co unt rate. To calibrate the SUMER data we perform the 



steps described in lMcIntosh et al.1 (|2007|) which includes corrections for flatfielding, as well as the 
absolute calibration. The last available flat-field is from 2001, which means that there are some 
uncertainties in correcting for detector artifacts. This is clearly visible in panel D of Fig. Q] which 
shows the SUMER Fe XII 1242 A intensity (calculated in the same way as described for EIS in the 
above). The uncertainty in flat-fielding is evidenced by the presence of bright and dark horizontal 
streaks in the raster. 

To coalign these rasters we perform the following steps. Since these rasters have very dif- 
ferent durations, they are subject to substantially different solar rotation effects with the features 
rastered by SUMER rotating by 30", with EIS only undergoing 10". To remove this effect, we 
recalculate the pointing information for each position of the EIS and SUMER rasters back to the 
same reference time t . In other words, we calculate which features on the Sun at to occupy the 
SUMER and EIS slits at the diverse range of positions and times of the rasters. Using this new 
pointing information, we then calculate which SUMER raster positions are nearest in space to 
the EIS raster positions, and form a SUMER raster that has the same spatial scale and number 
of positions as the EIS rasters. As a final step we perform a final coalignment between the new 
SUMER raster and the EIS rasters through visual comparison of the features in the eastern half of 
the overlapping parts of the field of view. We focus on the eastern half because the time difference 
between the SUMER and EIS exposures in that area of the field of view are the smallest, of order 
± 1 hour (see bottom panel of Fig. [3]), which increases the likelihood of viewing the same coronal 
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and TR structures. 



The resulting coalignment, shown in Fig. [H is quite satisfactory, especially in the eastern half 
of the raster. We use simultaneous TRACE and XRT data to identify the locations of TR moss. The 
high spatial resolution of the TRACE Fe ix/x 171 A image (formed at slightly lower temperatures 
of 1 MK, and taken around 17:45 UTC) shows the reticulated pattern that is typical for TR moss. 
This pattern is c aused by the presence of E UV absorbing chromospheri c features; typically jets o r 
dynamic fibrils JPe Pontieu et all l2003bfl . 120041 : lHansteen etaDbood : be Pontieu et alJl2007ah . 
We find several patches of moss around x = 20 — 40 and y = — 70 (MA), as well as x = 75 — 120 
and y = 60 — 90 (MB). Those footpoint locations fit very well with the location of the associated 
hotter loops in the two XRT images in Fig. [3l In addition, comparison of the Fe XII rasters with 
the XRT images allows us to identify coronal loops in the Fe XII rasters around x = 10 — 40 and 
y = 80 — 120 (A), as well as x = 60 — 120 and y = 30 — 60 (B), and some smaller loops around 
x = 10 - 25 and y = - 15 (C), x = 20 and y = 25 - 35 (D), x = 60 - 80 and y = 65 - 85 (E), 
and x = 40 - 65 and y = 60 - 75 (F). 

The XRT images (Fig. [3]) are taken at roughly the beginning and end of the SUMER raster 
and show the significant changes the western half of the active region undergoes during the 3 
hour raster. We draw attention to the formation of moss footpoints around x = 80 — 120 and 
y = 60 — 70 and associated loops just south of that region. These changes all occur after the EIS 
raster has ended. This explains the different appearance in panels A, B and D of Fig. Q] of the 
Fe XII emission in the southwestern part of the field of view. The eastern part of the active region 
shows much less variability with stable moss patterns. This is the region we will focus on in the 
following. 



3. Analysis of EIS/SUMER Observations 



We calculat e the density sensitivity of the Fe XII 186.88 and 195.1 A lines by using Chianti 
dDere et al.ll 19971) . The line ratio is sensitive to densities between log n e =8.5 and 10.5, as shown 
in panel F of Fig. \T\ We calculate the intensities as described in £j2[ and then determine the ratio 
of the 186.88 A and 195.1 A intensity. We use the curve in panel F to calculate the density shown 
in panel C of Fig. [TJ We find densities of order 10 9 cm -3 in the coronal loops we have identified 
earlier, with significantly higher densities of order 4 10 9 - 10 10 cm -3 in the moss regions. These 
values are very reasonable, and s i milar t o those found in ac tive region loops and moss regions 
reported by lFletcher & de Pontieu! (|1999i) : IWarren et al.l (|2008|) using CDS and EIS respectively. 



We can use the densities determined from this line ratio to calculate what the intensity ratio 
between the 195 A and 1242 A lines should be, given the density sensitivity of the latter ratio. This 
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density sensitivity is not as straightforward to use as the 186/195 ratio. First ly, there is some dis - 
cussion regardin g the atomic ph ysics underlying the density sensitivity, with lKeenan et al.l (I1990|) 
and CHIANTI (IDere et al.ll 1997H g i ving s lightly different curves of line ratio to density (see panel 
F of Fig. [T]), with the lKeenan etal.Nl990h results showing lower ratios by a factor of 0.73. In what 
follows we will us e both and compare t he results. Using the density derived from 186/195 A and 
the CHIANTI and iKeenan et al.1 (I1990|) curves, we now calculate for each position in the rasters 
the predicted Fe XII 195/1242 A intensity ratio (shown in panels H and I of Fig. [TJ). The values of 
this ratio vary between 8 and 12 for the CHIANTI curve, and 5 and 8 for the Keenan et al. dl990h 
curve. 

We now calculate the Fe XII 195/1242 A ratio directly from the observations using the coaligned 
rasters. The Fe XII 1242 A raster is quite noisy because of the low count rates, so we use a boxcar 
smoothing over 3 pixels to remove some of that noise when calculating the ratio between the EIS 
and SUMER rasters. The resulting ratio is shown in panel G of Fig. [T] with the same color table 
as panels H and I. The observed ratio varies between 5 and 15 and generally shows a much better 
correspondence with the levels predicted by the CHIANTI curve. However, it is immediately clear 
that there are also significant differences between the predicted and observed intensity ratios. The 
most striking difference is that the moss regions in the eastern half of the field of view system- 
atically show lower ratios (of order 4 to 6) than predicted (of order 8-10). This implies that the 
EUV emission with wavelengths below 912 A is weaker by a factor of order 2 in moss regions than 
would be expected. In contrast the loop structures we identified in the above show a much better 
correspondence with typically values around 10 (yellow) in both the observed and predicted maps, 
especially for loops A, C, D and E. 

Fig. [2] shows these differences in more detail by comparing histograms of the predicted and 
observed 195/1242 intensity ratio for coronal loops A and F (top panel), as well as moss region 
MA (bottom panel). As we can see, the loop regions show reasonable correspondence between 
the predicted and observed ratios, whereas the moss regions show a significant difference of order 
2. 

The loop regions often show higher observed values than predicted. This can be seen in Fig.[2l 
but also in the large loops B in Fig. [H We believe this is likely because the active region thermal 
structure in the western half of the field of view changed towards the end of the SUMER raster, so 
that the Fe XII 1242 A emission becomes much lower, perhaps because of a heating of the loops 
in that region towards the end of the SUMER raster (which is suggested by the XRT timeseries, 
Fig. [3]). In this scenario, the EIS raster, taken 3 hours earlier then still shows evidence of the cooler 
loops. 

Are the differences between the observed and predicted ratios caused by Lyman continuum 
absorption in the moss regions of 195 A emission? Our results hinge on the fact that the observed 
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spatial variation of the 195/1242 A ratio between moss and loops is larger than predicted by densi- 
ties determined from the 186/195 A ratio. We observe a generally much higher contrast and lower 
values in the map of observed ratios. Several mechanisms could potentially contribute to this mis- 
match between predicted and observed ratios. We argue in what follows that none of these can 
explain our observations and that Lyman continuum absorption is the most likely cause. 

A first effect is the temperatu re dependence of the 195/1 242 A line ratio. This temperature 



sensitivity was briefly discussed by lDel Zanna & Masonl (120051) and is illustrated in Fig. |4j We see 



from the top panel that for the range of densities we observe in this active region (log n e = 9 — 10) 
the 195/1242 A line ratio can vary from 5 to 3 for logT = 5.8 to from 15 to 10 for logT = 6.4. 
Both of these temperature values are on the extreme end of what we are most likely observing in 
these passbands, since the contribution functions of both lines are sharply peaked around logT = 
6.15. 

We believe that in the coronal loops this temperature dependence can explain some of the 
spatial variations of the 195/1242 A line ratio. In loops, we observe these variations to be larger 
than predicted. However, we argue that in the moss regions, the temperature dependence does not 
impact our result. We know from the 186/195 A ratio that the densities in the moss are of order 
log n e = 9.5 — 10. In principle our low observed values of the 195/1242 A ratio in the moss could 
thus be caused by a much lower temperature (log T = 5.8), since Fig. |4] shows ratios of order 3 
to 4. However, we know that moss is always associated with high pressure, hot and dense loops. 
So, if the upper transition region moss really had a temperature that is much lower than the peak 
of the Fe XII contribution function, then the high pressure, dense coronal loops that are associated 
with the upper TR moss would include plasma with temperatures at logT = 6.15. Since the loops 
associated with moss are high pressure and dense and at the temperature corresponding to the peak 
of the Fe XII contribution function, these loops would be very bright in the Fe XII rasters. However, 
we do not observe any such loops in or around the mossy regions. Our observations thus indicate 
that the temperature of the moss is not low enough to cause the low levels of 195/1242 A ratios 
that we observe in the mo ss. The lack of F e XII l oops connecting to moss has also been observed 



for other active regions by lAntiochos et al.l (120031) 



Another effect that could contribute to the mismatch between observed and predicted 195/1242 
A ratio in moss regions is the presence of blends. Any effects of blends can only impact our results 
and conclusion if the blends preferenti ally occur in loops or moss. We now show that the effect 



of blends does not change our results. lYoung et al.l (I2009|) have studied the density sensitivity of 



the 186.88 A and 195 A lines in detail. They find that there are several uncertainties in the de- 
termination of the intensities of the Fe XII 186 and 195 A lines that can potentially impact our 
analysis. There is a weak lin e at 186.98 A just redward of the Fe XII line, which is likely a Ni XI 



transition (| Young et alj|2009f) . In our analysis, we avoid contributions of this line by excluding the 
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wavelength range around the Ni XI line in our calculation of the total intensity of the F e XII 186.88 



A lin e. The Fe XII 186.88 A is also impacted by a weak blend of S XI at 186.84 A (|Young et al 



2009). Under equilibrium conditions, this li ne is formed at log T=6.3, slightly higher than the 
logT = 6.15 temperature of the Fe XII line. lYoung et al. (i2009h find that the S XI intensity is 



typically 2-5% of the Fe XII 186.88 A line, which is smaller than the photon noise on the observed 
Fe XII 186.88 A intensities. Given the density sensitivity of the 186/195 A ratio (Fig. [D, such a 
small difference cannot lead to a significant reduction in the discrepancy between predicted and 
observed 195/1242 A intensity ratio in the moss regions. To remove the factor of 2 discrepancy 
in 195/1242 A ratio would require a change in the 195/186 ratio by a similar factor. Clearly the S 
XI blend is too weak to impact our results. In addition, even if the blend was much stronger than 
previously observed, it would more likely impact the hotter moss regions, because the S XI line 
is formed at slighly higher temperatures. That would mean that in the moss the Fe XII 186.88 A 
intensity is overestimated by our analysis, which would lead to a predicted 195/1242 A intensity 
ratio with even less difference between moss and loops than is observed. 



Young et al] (120091) shows that the Fe XII 195.12 A line similarly has some weak lines in its 
immediate vicinity. We again avoid contributions of thes lines by excluding the wavelength range 
of these weaker lines in our calculation of the total intensity of the Fe XII line. This line is also 
blended by another Fe XII line, which is at levels of less than 10% for the densities we find here 
(< 10 10 cm -3 ). Given the weakness of the blend, and the fact that the brightness of the blend 
should not show any difference between moss and loops (since it is also Fe XII), this blend cannot 
change our conclusions. 



Young et al.l (|2009|) have compared densities from different line pairs and found some dis- 
crepancies. They suggested that the density determination from the Fe XII 186.88 and 195 A pair 
may lead to overestimates of the density by 0.2-0.4 dex, especially in higher density regions. Since 
our observations show that the moss regions have the highest densities, this currently unexplained 
effect (most likely due to uncertainties in the atomic data) could lead us to overestimate the moss 
densities. However, overestimates of the density in moss would lead to even higher values for the 
predicted ratio of Fe XII 195 and 1242 A intensity. In summary, this issue would not remove, but 
rather exacerbate the discrepancy between observed and predicted ratios. 

We have also investigated the effects of the main sources of noise on the EIS spectra. We 
calculated the errors on our intensities using the eis_prep software that is part of the EIS solarsoft 
tree, and found that photon (Poisson) noise dominates the relative errors. In the moss and loop 
regions we study, the spectra are high signal to noise, of order 20 for the 186.88 A line and 30- 
40 for the 195 A line. This means that the relative errors o on the intensities are small, of order 
5 %. If we assume a worst case scenario in which, e.g., the 186.88 A intensity is lowered by 
cr 186 and the 195 A intensity is increased by a 195 , the predicted 195/1242 ratio map is shifted to 
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higher values. However, such a shift occurs in equal measure for moss and loop regions, so that 
the predicted spatial variation of the 195/1242 A ratio cannot be reconciled with the observed 
ratio. In other words, while these errors can introduce noise on the predicted ratio, they cannot 
explain the consistently lower values of the 195/1242 ratio that we observe in the moss regions. A 
similar argument can be made for the absolute calibration errors between SUMER and EIS. Such 
errors will not depend on the nature of solar region observed (especially since those are of similar 
brightness), and cannot explain the spatial variation we observe. 

We conclude that the most likely explanation for our EIS/SUMER observations is that EUV 
emission (at 195 A) in upper TR moss is reduced in brightness by a factor of order 2 because of 
absorption from neutral hydrogen and helium, and singly ionized helium in the chromospheric jets 
that permeate the moss. 



4. Analysis of STEREO Observations 

The region we study with EIS and SUMER was observed by Hinode and SOHO close to disk 
center. However, it was also observed by both STEREO spacecraft (SC). The multiple vantage 
points these spacecraft provide allow us to also investigate the effect of the angle 9 between the line- 
of-sight vector and the local vertical on the chromospheric absorption of TR EUV emission. We use 
simultaneous STEREO A and B images taken in the Fe XII 195 A passband of SECCHI on 14-Nov- 
2007 during the EIS and SUMER rasters (i.e., between 16 and 18 UTC). The images were corrected 
for dark current using the software provided in the STEREO solarsoft package. To account for the 
different distance to the Sun of both spacecraft, we correct the intensities observed by STEREO 
A by multiplying the SC A intensities with a factor of (Rsuua/ Rsuub) 2 = (994/920) 2 , where 
RsuriA and Rsun B is the solar radius in arcsec as seen from SC A & B respectively. 

We study the brightness of three regions of interest in this active region as observed with SC 
A and B. Fig. |5] shows a SC B image of the two moss regions we selected as well as a loop-like 
structure. All of these regions also occur in the EIS and SUMER rasters of Fig. \T\ The region is 
seen at disk center by STEREO A, and at an angle of 40 degrees between line-of-sight and local 
vertical by STEREO B. We determine a histogram of intensities seen by SC A and by SC B for 
the three regions of interest. We find that the average intensity of the loop is the same for SC A 
and B (to within 5%). The observed discrepancy is within the range expected from the absolute 
calibration error, which is estimated to be of order less than 10% [± In contrast, the histogram of 
intensities for both moss regions (Fig. [6]) shows significant differences, with the emission seen by 



1 Wuelser, 2008: private communication 
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SC B reduced to 75% of that seen by SC A. Moss region A shows an average intensity of 340 
(normalized units) with a standard deviation of 130 when seen by SC A, whereas it is weaker at 
240±93 for SC B. Similarly, moss region B is at 310 ± 100 for SC A, and 260 ± 60 for SC B. This 
is confirmed by the histogram plots, which show a marked increase of very bright features seen 
with SC A (Fig. [6]). 

We find that the coronal loop intensities do not change, whereas the moss region intensities 
decrease significantly when viewed from the side (at a viewing angle of 40 degrees). These results 
indicate that there is a significant increase of absorption by neutral hydrogen or helium in moss 
regions when those are viewed from the side. The fact that the loop intensities do not change 
significantly is to be expected since such loops occur high in the corona, where there is little 
neutral hydrogen or helium. In contrast, the mix of cool jets and hot plasma in moss regions 
implies that an observer who views the moss from the side (at a non-negligible viewing angle 
from the local vertical) will see a larger surface area of chromospheric plasma, which will lead to 
increased absorption. The amount of absorption seems to roughly scale, perhaps coincidentially, 
as cos 9, with 9 = 40 degrees for our observations. 



5. Analysis of Simulations 

The presence of significant chromospheric absorption of TR EUV emission with wavelengths 
shorte r than 912 A is confi rmed by advanced radiative MHD simulations. In the Oslo Stagger 



Code (IHansteen et al.ll2007|) models of the solar atmosphere stretching from the convection zone 
to corona are constructed. These models include convection, non-grey, non-LTE radiative losses 
in the photosphere, chromosphere, and corona, as well as conduction along the magnetic field. 

An initially potential magnetic field is inserted into a fully convective model of dimensions 
16 x 8 x 15.5 Mm 3 of which 1.5 Mm is below the photosphere and 14 Mm above. The grid size 
in the horizontal direction is 65 km, with a smaller grid size of 32 km in the vertical direction. 
Convective motions and photospheric granular dynamics stress the magnetic field resulting in an 
upwardly propagating Poynting flux that is dissipated in the outer solar layers. After approximately 
20 minutes solar time, and depending on the initial magnetic field strength and topology, the dis- 
sipated energy is sufficient to raise and maintain the temperature in the upper layers to coronal 
values. Thus a simulated atmospheric structure emerges, with photosphere, chromosphere, and 
corona much as one believes exists on the sun. In the present context we point out that the chro- 
mosphere — corona interface, the transition region, in this model is extremely corrugated ranging 
in height from some 1.5 Mm to 5 — 6 Mm above the photosphere. Thus at any given time we 
find plasma at temperatures ranging from 5, 000 K to 1 MK in the same height range. As a conse- 
quence cool plasma in sufficient quantities to absorb short wavelength EUV radiation is found at 
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these heights. 

To quantify this we calculate the intensity in both the Fe XII 1242 A and 195 A lines from a 
given snapshot. This is accomplished by integrating the optically thin contribution function 

dl = A c \n 2 c g{T)e~ T dz 

along a rays parallel to the "line of sight", where A e \ is the element abundance relative hydrogen, 
n e is the electron density, r is the optical depth of the absorbing opacity, and g(T) is a function 
that incorporates the ionization state of the radiating material and the collisional excitation rates: 

g(T) = O.SShvnjAjifW/ne. 

Here hv is the energy of the transition, nj is the population of the upper level of the radiating 
ion, Ajj the Einstein coefficient for the transition, and f(T) the ionization state of the radiating 
ion (|Dere et al.lll997l : LLandi et al J 120061) . The optical depth r is computed assuming that all ab- 
sorption comes from ne utral hydrogen and neutral and singly ionized helium such that following 
Anzer & Heinzell (|2005l) we write 



T = Or I-Wr I + He I^He II + He llNue II 

where i\T H i is the total column density of neutral hydrogen along the line of sight, and N^ e i 
and iVne n are those of neutral and singly ionized helium. The o-values are the respective pho- 
toinization cross sections which depend only on wavelength . In computing t he col umn densities of 
hydrogen and helium we assume the ionization derived by iMazzotta et al. (1 19981) as found in the 
CHIANTI package. 

In figure [7] we plot the emergent intensities of the Fe XII 1242 A and 195 A lines as seen 
from the side, e.g. as seen at the limb. Notice the significant differences between the intensities, 
both in the higher-lying loops (because of density sensitivity), and at the lower regions near the 
footpoints of the loops (z < 3 Mm) because of the chromospheric absorption of EUV photons that 
have wavelengths below 912 A. The presence of almost complete extinction of the 195 A intensity 
in the lower 1-3 Mm of the coronal loops is caused by absorption from upper chromospheric 
extrusions into coronal territory. 



6. Summary 

We find evidence for significant absorption of EUV emission (with wavelengths below 912 
A) in the transition region of coronal loops because of the presence at similar heights of dense, 
chromospheric gas. We use the Fe XII 186.88, 195 and 1242 A lines from quasi-simultaneous 
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Hinode/EIS and SOHO/SUMER spectra of an active region and show that the transition region 
emission in moss regions is reduced by a factor of order 2 because of absorption by neutral hydro- 
gen and helium and singly ionized helium that occurs at similar heights in chromospheric jets. We 
use STEREO/SECCHI observations of the same active region using both spacecraft and find that 
this absorption is further increased when moss regions are viewed closer to the limb. These obser- 
vational results are confirmed by advanced 3D numerical simulations which show similar levels of 
absorption because of co-mingling of cool, chromospheric plasma and hot, upper transition region 
plasma at similar heights. 

The determination of this type of absorption significantly reduces or removes an uncertainty 
that has plagued comparisons of EUV images with numerical models of coronal loops, which have 
had trouble reconciling the observed and predicted EUV brightness at the footpoints of loops. Cur- 
rent models typically show a mismatch of one order of magnitude. A variety of attempts at reduc- 
ing the discrepancy include height (and magnetic field) dependent filling factors, time-dependent 
heating models, etc... By determining what the real, unobstructed emission in the upper transition 
region of coronal loops is, our results reduce (but do not totally remove) the current discrepancy 
and will guide modelers towards better constraints on their free parameters and assumptions. The 
STEREO observations and numerical results suggest that full Sun models need to also consider the 
center-to-limb variation of the absorption. 

The significant absorption of EUV emission in the transition region by cool plasma in chro- 
mospheric jets indicates that coronal loop modellers should be aware of and take into account the 
properties of the absorbing jets. For example, is it possible t hat the observed inv erse correlation 



between "filling factor" and coronal loop pressure reported by [Warren et al.l (|2008|) is related to the 
dependence on coronal pressure of the height range over which the transition region occurs in com- 
bination with the presence of absorbing jets? Given the significant absorption by chromospheric 
jets, more attention is also warranted for the effects of chromospheric absor ption on the detailed 



shape of line profiles of coronal emission lines (see, e.g. lSpadaro et al.l l 1996 s ) 



The total amount of absorption we find is of order at least 2 to 3 (for viewing angles of 
less than 40 degrees from disk center). Further observations with STEREO/SECCHI for a wider 
variety of viewing angles have the potential of revealing the true center to limb variation of moss 
brightness, which can be compared to that of numerical simulations, and may reveal details about 
the filling factor and orientation of chromospheric features in the moss phenomenon. In addition, 
simultaneous observations with EIS and SUMER of a larger statistical sample of active regions 
will help determine the variability of this absorption across active regions of varying activity. 
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Fig. 1. — Top row shows a spectroheliogram in Fe XII 195.1 A and 186.88 A taken on 14-Nov- 
2007 from 16:44 to 17:50 UTC with EIS, as well as the logarithm of the density as determined 
from the line ratio of the 186.88 and 195.1 A lines. The middle row shows the a slightly spatially 
smoothed SUMER raster taken from 17: 12 through 20:33 UT and a TRACE Fe ix/x 171 A image. 
The middle right panel shows the line ratio to density functions we used to determine the density 
for both line pairs. The Fe XII 186/195 A line ratio is multiplied by 10 for illustration purposes. 
The bottom row shows the observed Fe XII 195/1242 A line ratio, as well as the line ratio that 
is predicted if we use the density de t ermin ation of the 186/195 line ratio, for two atomic models: 



CHIANTI and that of Keenan et al. 



1990J). The overplotted contours correspond to a brightness 



threshold in the EIS Fe XII 195.1 Araster. See text for details on the interpretation. 
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Fig. 2. — Comparison of histograms of observed (full lines) and predicted (dotted lines) values of 
the 195.1 A to 1242 A intensity ratio in coronal loops A and F (top row), and moss region MA 
(bottom row). See § [2]for the definition and location of loop and moss regions. 
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Fig. 3. — Top two images show XRT images (Al po i y ) taken at the beginning and end of the SUMER 
raster. The overplotted contours again correspond to the same brightness threshold in the EIS 
Fe XII 195 A raster shown in Figure [TJ Since the XRT filter is sensitive to hotter plasma than the 
Fe XII rasters, we see much more emission from the hot loops connecting the TR moss footpoints. 
In addition, we see that the loop structure in the middle right of the image drastically changes in 
morphology with new moss footpoints appearing around x = 80 — 120, y = 60. The lower panel 
shows, in hours, the time difference between the SUMER and EIS exposures (shown in Fig. [[]) at 
each coaligned raster location. The time differences are caused by the different durations of the 
SUMER and EIS rasters, and the fact that they raster in E-W and W-E directions respectively. The 
left part of the rasters is taken within 1 hour of each other, whereas the right part (x > 70) shows 
exposures that are taken 2 or more hours apart. See text for details. 
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Fig. 4. — The density sensitive intensity ratio of the Fe XII 195 A and 1242 A lines also varies 
depending on the temperature of the plasma (form the CHIANTI atomic database). The left plot 
shows the density sensitivity for a range of temperatures with logT from 5.8 through 6.4. The 
temperature associated with the maximum formation is logT = 6.15, with the contribution func- 
tion dropping off steeply towards lower and higher temperatures. The right plot shows the same 
dependency, but now the ratio is plotted as a function of temperature for a range of densities. The 
densities in the active region studied vary between logn e = 9 — 10. 
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STEREO B: 195 Channel 
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Fig. 5. — Fe XII A image of the same active region studied with EIS and SUMER, as seen with 
STEREO B. Marked are three regions of interest (one loop region, and two moss regions) for which 
the intensities seen with STEREO A and B are compared (see Fig. [6]). STEREO A observed this 
region at disk center, whereas STEREO B observed the region at a viewing angle of 40 degrees 
(with the local vertical). 
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Fig. 6. — Histograms of the intensity in the two moss regions (marked as A and B in Fig. [5]) for 
Fe XII 195 A observed with STEREO A (full line) and STEREO B (dashed line). Both regions 
are reduced in intensity by a factor of roughly 0.75 when viewed with STEREO B. This reduction 
in intensity is caused by Lyman continuum absorption of EUV emission by chromospheric plasma 
that occurs at similar heights as the EUV emission. Such absorption is increased when the region 
is viewed from the side, since the observable surface area of chromospheric plasma is increased. 
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Fig. 7. — Side views of synthetic Fe XII intensities from our 3D radiative MHD simulation. The 
intensity shown in each panel is summed along a ray parallel to the y-axis (perpendicular to the 
plane of the figure) using the density and temperature in the simulated box, and CHIANTI to 
calculate the emergent intensity. The top panel shows Fe XII 1242 A emission, whereas the bottom 
panel shows Fe XII 195 A emission. Notice the significant differences between the intensities, both 
in the higher-lying loops (because of density sensitivity), and at the lower regions or footpoints of 
the loops (z < 3 Mm) because of the absorption of EUV photons that have wavelengths below 912 
A. See text for details. 



